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Within the framework of macroscopic quantum electrodynamics, the resonant van der Waals po- 
tential experienced by an excited two-level atom near a planar magneto-electric two-layer system 
consisting of a slab of left-handed material and a perfect mirror is studied. It is shown that disre- 
garding of material absorption leads to unphysical results, with divergent values for the potential 
away from the surface. Under appropriate conditions, the setup is found to feature a barrier near the 
surface which can be employed to levitate particles or used as a trapping or cooling mechanism. Fi- 
nally, the problem of spontaneous decay [J. Kastel and M. Fleischhauer, Phys. Rev. A 68, 011804(R) 
(2005)] is revisited. Disregarding of absorption is shown to drastically falsify the dependence on the 
atomic position of the decay rate. 

PACS numbers: 12.20.-m, 42.50.Wk, 34.35. +a, 42.50.Nn 



I. INTRODUCTION 

The first most comprehensive account of the proper- 
ties of left-handed materials (LHMs), i.e., materials si- 
multaneously exhibiting negative real parts of the elec- 
tric permittivity and the magnetic permeability in a cer- 
tain frequency interval, was given by Veselago [l|. As 
he pointed out, a negative refractive index may be at- 
tributed to such materials, leading to a number of un- 
usual optical phenomena. While natural materials with 
a negative real part of the permittivity such as plasmas Q 
or a negative real part of the permeability such as some 
ferrites 3] exist, LHMs have not yet been found in na- 
ture. During the last decade, significant progress in the 
fabrication of LHMs in laboratories has been achieved. 
Metamaterials based on periodic arrays of split ring res- 
onators and a grid of metallic rods have been reported 
to exhibit left-handedness in the microwave region 0, Q . 
More recently, a two-dimensional material consisting of 
a metal-insulator-metal waveguide has been experimen- 
tally demonstrated to show a negative refractive index in 
the optical region @ . The effective permittivity and per- 
meability of such metamaterials, valid on length scales 
which are sufficiently large in comparison to the elemen- 
tary building blocks of the material, are determined ei- 
ther theoretically 0, W\ or by means of reflection experi- 
ments 0. 

As already indicated in Ref. 1], a planar LHM slab 
with a negative refractive index has the interesting prop- 
erty of being able to focus light — an effect which suggests 
the possibility of realizing a so-called superlens for sub- 
wavelength imaging with a resolution well beyond the 



diffraction limit JJ. The superlens concept has been a 
subject of intense discussion [13, [lH, and limiting fac- 
tors such as the finite dimension of the lens [I4I or the 
influence of absorption have been studied [ij] . 

In view of the unusual properties of LHMs, it has also 
been of interest to study the modifications of dispersion 
forces in the presence of magneto-electric materials. In 
particular, the question of whether a left-handedness may 
lead to a repulsive Casimir force between two magneto- 
electric plates has been addressed [Tj]. From the study 
of van der Waals (vdW) forces on single [la] and be- 
tween two ground-state atoms '16' in the presence of 
a planar system exhibiting magneto-electric properties, 
it has been shown that the attractive/repulsive behav- 
ior of the potential is influenced by the strength of the 
magnetic properties rather than by the left-handedness. 
This may be attributed to the fact that in both cases 
the atoms are in the ground state, so the vdW forces 
represent off-resonant interactions that involve integrals 
over the whole frequency range. Hence, frequency in- 
tervals where no left-handedness occurs will inevitably 
also contribute to the force and counteract the effect of 
the left-handedness; this becomes evident after express- 
ing the force in terms of an integral over the imaginary 
frequency axis where both the permittivity and the per- 
meability are always real and positive. Similar arguments 
also apply to the abovementioned Casimir force between 
two magneto-electric plates, where a repulsive behavior 
has also been found to primarily require strong magnetic 
rather than left-handed properties [l4| . 

Two strategies might be envisaged to enhance the in- 
fluence of left-handed properties on dispersion forces. 



The first one would be to consider a material exhibiting 
left-handedness over a sufficiently large (real) frequency 
rangej_as recently discussed in the context of the Casimir 
force IIT"!; however, such an requirement is, in general, in 
conflict with the causality requirements met by all phys- 
ical media. Secondly, one may consider excited systems 
like amplifying media as also addressed in Ref . [17] . The 
fundamental difference between absorbing and amplify- 
ing media requires, however, a careful reconsideration 
of the underlying quantum electrodynamics [l8|. The 
situation is more transparent in the case of an excited 
atom where the necessary theory of vdW forces is read- 
ily available [l^. Excited atoms give rise to resonant 
force components which directly depend on the frequen- 
cies of possible downward transitions. If a material ex- 
hibits left-handed properties at one such frequency, one 
may expect these to influence the dispersion forces in a 
more significant way. 

The spontaneous emission of an excited atom situated 
in vacuum in the vicinity of a perfect mirror combined 
with a lossless LHM slab is studied in Refs. [20, l2]| . 
where it is noticed that a slab of unity negative refrac- 
tion (n = — 1) has the peculiar feature of making the space 
between atom and mirror appear of zero optical length 
when the atom-slab distance is equal to the thickness 
of the slab. For such positions, called focal points, com- 
plete suppression of the spontaneous emission for a dipole 
moment parallel to the surface, and enhancement of the 
spontaneous emission rate by a factor of two for a dipole 
moment perpendicular to the surface is reported. More 
recently, it has been pointed out that in reality inhibition 
of spontaneous decay can be weakened due to nonradia- 
tive decay at short distances and due to radiative decay 
at large distances 22|. 

Motivated by Refs. [13, HH, we shall study the vdW 
interaction of an excited two-level atom with the quan- 
tized electromagnetic field in the presence of a perfect 
mirror combined with a realistic LHM slab. To conform 
with causality and to avoid unphysical consequences, ab- 
sorption in the LHM is taken into account from the very 
beginning. The paper is organized as follows. In Sec. m 
basic formulas for the vdW potential and the decay rate, 
and for the Green tensor of a planar magneto-electric 
multilayer structure are given. Section IIIII is devoted to 
analytical and numerical studies of the vdW potential 
experienced by an excited two-level atom in front of a 
system consisting of a LHM slab combined with a per- 
fect mirror. Special emphasis is put on the effects of 
material absorption. In Sec. IIVI the atomic spontaneous 
decay rate is discussed and a comparison with earlier re- 
sults is made. Some concluding remarks and a summary 
are given in Sec. IVl 



eigenstate |n), transition frequencies ujnk, electric-dipole 
transition matrix elements d„fe, polarizability Q!(tt')] in 
the presence of an arbitrary magneto-electric medium 
of complex permittivity e(r,Lu) and permeability fi(r,Lu). 
The medium affects the classical Green tensor G(r, r', oj) 
of the electromagnetic field, which satisfies the differen- 
tial equation 
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G(r,r',a;)=<5(r-r') 
(1) 



/i(r,w) 
together with the boundary condition at infinity 

G(r,r',w)^0 for |r-r'|^oo. (2) 



The potential can be decomposed into a resonant and an 
off- resonant part (see, e.g., Ref. [l9l|). 



C/„(r^) = C/;(r^)-fC/-(rA), 
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(0: unit step function) and 
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respectively. Here, G*^^-'(r,r',Ci;) refers to the scattering 
part of the Green tensor. For an excited atom, the res- 
onant part of the potential often dominates and the off- 
resonant part can be neglected, as we will do throughout 
this paper. We shall also be interested in the rate of the 
spontaneous decay \n) ~* \k) of the atom, which reads as 
(see, e.g., Ref. Q) 



F = Q{uJnk) 
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d„fc • ImG(rA,rA,Wnfe) • d„fc. (6) 



In what follows, we restrict ourselves to a two-level 
atom with upper level |1) and lower level |0). Consider 
the atom to be placed in free space in front of a meta- 
material slab of thickness d, permittivity e{io) and perme- 
ability /i(ti') which is bounded by a perfectly conducting 
mirror on the back, as sketched in Fig.[TJ The coordinate 
system is chosen such that the z-axis is perpendicular to 
the slab and its origin coincides with the slab-vacuum 
interface. The scattering part of the Green tensor at the 
relevant atomic transition frequency and equal positions 
r = r' = r^ in the vacuum region is given by [15J 



II. BASIC EQUATIONS 

Consider the vdW potential of a neutral, unpolar- 
ized and nonmagnetic atomic system [position r^, energy 



G^'^^{rA,rA,UJnk) = G'^^'' {ZA,ZA, i^io) 
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In Eq. ([7]), p (s) denotes p- (s-)polarization. The reflec- 
tion coefficients are given by 
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with 
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where the wave vectors in the z direction are given by 



(10) 



/3=\/=^-9'' Pi = Jkl-q^. 



(11) 



with A:i, which always appears in the form of fc^ in the 
Green tensor, Eq. ([7]), being given according to 




FIG. 1: Atom in a front of a slab of an (absorptive) nieta- 
material of thickness d wliich is bounded by a perfect mirror 
on the back. 



FIG. 1; Atom in a front of a slab of an (absorptive) meta- 
material of thickness d which is bounded by a perfect mirror 
on the back. 
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(12) 



The s- and p-polarization unit vectors are defined as 

e^^GgXe^, e^ =: (ge^=F/3eg) (13) 

^10 

(e, — (\/q, q = |q|). The square root in Eq. (fTT|) has 
to be chosen such that it obeys the physical requirement 
Im/3i > for a passive medium, with this choice being 
unique for purely absorbing media (Ime,Im/i > 0) [8|. 

(Reelm/z + 



The imaginary part of Pf reads Im/3j 
ImeRe/i)(jJ2o/c^- It is positive for right-handed materi- 
als, which together with the condition Im/3i > implies 
that /3i has to be in the first quadrant of the complex 
plane. On the other hand, Im f3f is negative for left- 
handed materials. This, combined with the condition 
Im/?i > 0, leads to that /3i lies in the second quadrant 
of the complex plane (see also Ref. |24|V 

Introduction of polar coordinates in the g-plane, i.e., 
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and, with e^ = (=F/3ccos0/ci;io, =F/3csin(/)/a;io, ^c/wio), 

(3'^ cos^ 4> —P'^sm(l)cos(l)—q(3cos(f)\ 
'fi'^ sin^ q 
qP sin (j) 
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(15) 
With d^q — qdqdcj), the angular integration in Eq. ^ 
can be performed on using 
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resulting in 
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From Eq. ^TE\i it can be seen that an atom with a dipole 
moment perpendicular to the surface is coupled to the 
p-polarized waves only while an atom with a dipole mo- 
ment parallel to the surface is coupled to both p- and 
s-polarized waves. It is instructive to decompose the in- 
tegral in Eq. (fT5)) into two parts, 
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The first integral, which contains an oscillating factor, 
results from propagating waves, whereas the second one, 
which contains an exponentially decaying factor, results 
from evanescent waves. The first is the more important 
one away from the surfaces while the second dominates 
near the surfaces. 



III. VDW POTENTIAL 

A. Ideal (nonabsorbing) LHM 

Let us assume hypothetically an absolutely nonabsorb- 
ing LHM having e = /i = — 1 . In accordance with the re- 
marks below Eq. p3)) . we thus set 



/3i 



-/? for q < Lo/c, 
j3 for q > Lu/c, 



and the reflection coefficients 
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Note that this result does not depend to the sign of the 
square root chosen for (3i, since for the ideal LHM dis- 
cussed here, the reflection coefficients are invariant under 
a change (3i — > — /3i. Substitution of these reflection co- 
efflcients in Eq. (|18p for the Green tensor leads to 
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(22) 



Equation ([22|) reveals that for za > d, the Green tensor 
exactly coincides with that of a conflguration in which 
a perfectly conducting mirror is placed at z = d. As 
can be seen from an image dipole construction, this is 
due to the perfect negative refraction taking place at the 
vacuum-LHM interface: The image of an electric dipole 
ai ZA > d is situated at z^ = d — {za — d), just as it 
would be if a perfectly conducting mirror were placed at 
z — d (see Fig. [2]). An important difference to the case 
of a mirror placed at z = d is the fact that in the setup 
considered here the atom can be placed in the region 
Q < ZA < d. However, it can be seen that, whereas the 
flrst integral in Eq. (|22p , which is complex, is always well- 
behaved, the second (purely real) integral tends to minus 
infinity in this region. Thus we are left with an infinite 
potential for < Zyi < d, despite the absence of any 
physical surface at d. We clearly have a situation where 
an unphysical assumption (of absolute zero absorption) 
gives rise to unphysical results. 

After calculating the two integrals, we obtain for 
za> d 



G^xl {za, za, ujw) = G,[,^J {zA,ZA,uJw) 
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FIG. 2: Image-dipole construction for the setup depicted in 
Fig. [T] The dashed line marks the position of a perfect mirror 
that would generate the same image. 
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[z = 2ujio{za — d)/c]. Note that the off-diagonal com- 
ponents of the Green tensor are equal to zero. For 
ZA > d, the vdW potential then reads, in accordance with 
Eqs. ©-(O and after neglecting the off-resonant part. 
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[^w = {{dw)^,{dio)y,Qi), djL„ = (0,0,(dio)^)], where 



ReG^ = -^ [cos(z) -K5sin(i) - 52cos(5)] (26) 



and 
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For Q < Za < d, the potential is, as already mentioned, 
divergent due to the unphysical assumption of vanishing 
absorption. Hence the question may arise whether at 
least in the range za> d the above-given result can be 
regarded as being correct in the limit of sufficiently small 
absorption. An answer to this question can obviously 
only be given by taking material absorption into account 
from the very beginning, as shall be done in the next 
section. 



B. Real (absorbing) LHM 

To allow for material absorption, we return to Eq. (jlSp 
for the scattering part of the Green tensor and set therein 
e(wio) = — 1 -I- irje, /i(a'io) = — 1 -I- irim, where for sim- 
plicity we assume rj = rje ^ rjm- Due to the positive 
imaginary parts of e and /i, divergent integrals of the 
type of the second integral in Eq. ((22|) can never occur. 
Note that nonvanishing absorption helps to regularize the 
behavior of the integrals, just like it does in the superlens 
geometry [3|. 



In Fig. [3l the resulting (resonant) vdW potential is 
plotted versus the distance from the atom to the LHM 
slab, for the two cases of parallel and perpendicular align- 
ment of the atomic dipole moment with respect to the 
slab and for different values of absorption. It is seen that 
the potential features an attractive behavior in the non- 
retarded regime, which is governed by an inverse power 
law, while in the retarded regime an oscillating behavior 
with alternating sign of the potential occurs. The non- 
retarded regime is dominantly determined by evanescent 
waves [first integral in Eq. (|19p ] while propagating waves 
are mainly responsible for the retarded regime [second 
integral in Eq. ^9]) ]. The figure further reveals that in 
the limit of vanishing absorption when 77 tends to zero, 
then the potential approaches the one found, for za > d, 
in the case of absolutely nonabsorbing LHM [Eq. pS)) 
together with Eqs. (|26p and ((27|) ]. It is clearly seen that 
in practice unavoidable absorption always prevents the 
system to feature an infinitely negative potential in the 
interval < za < d, as would be expected from the 
study of a nonabsorbing LHM. We do, however, note 
that for very small absorption, the potential starts to be- 
come strongly negative around the region za ~ d, where 
this effect is more noticeable on the case of perpendicular 
atomic dipole moment. 

A very striking feature of the potential is completely 
missing if absorption is disregarded: For a transition 
dipole moment parallel to the surface and medium ab- 
sorption, a barrier appears in front of the slab at dis- 
tances ZAWio/c < 1 [Fig. [3i;a), 10-4 < ?? < 10~^]. The 
fact that barriers occur only for transition dipole mo- 
ments parallel to the surface but not for those perpen- 
dicular to the surface suggests that s-polarized waves, 
which are coupled to the first but not the latter, play an 
important role in their formation. 

For sufficiently weak absorption [rj = 10^^ in Fig. ^a.) 
and 77 = 10-3, IQ-'', 10"^ in Fig. M.^)], an attractive po- 
tential starts to appear at distances of a few wavelengths, 
somewhat away from the surface. Atoms located within 
this range will get adsorbed to the surface. This behav- 
ior is more pronounced for a transition dipole moment 
perpendicular to the surface than for a transition dipole 
moment parallel to the surface — a fact that can be re- 
lated to the coupling respective no coupling of these two 
dipole moment orientations to s-polarized waves. By set- 
ting the denominators in Eqs. ^ and © equal to zero, 
one finds that the poles of the generalized reflection co- 
efficients for p- and s-polarized waves, respectively, are 
solutions to the equations 
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FIG. 3: Resonant VdW potential experienced by an excited 
two-level atom in the setup sketched in Fig. [1] for d = 5c/a;io, 
£(ljio) = /i(tJio) = —1 + 177, and (a) dipole moments parallel 
to the surface and (b) dipole moments perpendicular to the 
surface. Curves 1-5 refer to various degrees of absorption 
T] = 10-\ 10-^ 10"'*, 10"^ and 0, respectively. 



These are the dispersion relations for the surface plas- 
mon polaritons [25| , which can be strongly excited when 
the atom is close to the surface. The real parts of the 
poles lie in the regions covered by the second integral in 



Eq. ((T9| . The poles of the s- (p-)polarized waves are in 
the upper (lower) half of the complex plane, where for 
smaller absorption they are closer to the real axis. 

The influence on the vdW potential of the slab thick- 
ness is illustrated in Fig. U) Note that zero slab thick- 
ness means an atom in front of a perfect mirror. As the 
slab thickness increases, a potential barrier arises and 
grows in height for dipole moments parallel to the sur- 
face. However, at some threshold value of d, the height 
of the barrier starts to be reduced, and eventually the 
barrier disappears when the slab is too thick. This can 
be explained as resulting from the increasing effects of 
material absorption. It is seen from the figures that the 
distance from the surface at which a relatively strong at- 
tractive potential can occur increases with the slab thick- 
ness. This is more pronounced for perpendicular dipole 
moments rather than for parallel ones — a consequence of 
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FIG. 4: Resonant vdW potential experienced by an excited 
atom in the setup sketched in Fig. [T]for e(a;io) = /i(wio) = 
— 1 + ilO~^, and (a) dipole moments parallel to the surface 
and (b) dipole moments perpendicular to the surface. Curves 
1-5 refer to various thicknesses of the slab du\o/c = 0, 3, 5, 
7 and 10, respectively. 



the fact that, as already mentioned, the first is coupled 
more efficiently to the evanescent waves. 



C. Near-surface limit 

Further insight into how the appearance of the barrier 
depends on the amount of absorption, the orientation 
of the atomic dipole moment as well as the thickness of 
the LHM slab can be gained by examining the vdW po- 
tential in the near-surface limit. Near the surface, the 
evanescent waves dominate the potential [the second in- 
tegral in Eq. (fT9|) ]. Let us consider the (resonant part of 
the) potential as given by Eq. ^ together with Eqs. (fT5|) 
and (fTO|) for very short distances, za^io/c <C 1, in more 
detail. The main contribution to the g-integral is from 



values q » ujw/c and q » ^J\£^,\LUlQ/c [e = ^(wio), 
/i = /^(wio)], since the exponential factor [in the sec- 
ond integral in Eq. (fT9|) ] effectively limits the integration 
interval to values q"^ za- Hence, with 



(3 ~ iq, /3i ~ iq 
the potential in the short-distance regime reads 
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where 



RerL = [(|e|2-l)(l + e-4^'^) 
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then the first terms in Eqs. ([5^ and ([55]) approximately 
vanish, thus 
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where the opposite signs imply that the two polarizations 
give competing contributions to the potential. Namely, 
the p-polarized waves give rise to attractive contributions 
to the potential while the s-polarized waves lead to repul- 
sive ones. Very close to the surface, due to the presence of 
the q^ factor [see Eq. ((3T|) ] , the contribution to the poten- 
tial of the p-polarized waves is proportional to ^j z\ while 
the contribution of the s-polarized waves is proportional 
to \l ZA- The contribution of the p-polarized waves hence 
dominates, resulting in an attractive potential (also see 
Figs. [3] and [4|). At some distance from the surface, the 
contribution of the s-polarizcd waves can dominate under 
appropriate conditions, which then leads to the appear- 
ance of a potential barrier. This also explains the absence 
of the barrier in the case where the dipole moment is per- 
pendicular to the surface. 

Inspection of Eqs. (|32p and 
is positive for all |£| > 1 while Rer2_ is negative for all 



reveals that Rer2_ 



|/i| > 1. Since the height of the potential barrier is de- 
termined by the magnitude of lRer|_|, lefthandedness is 
not necessarily required to realize a potential barrier in 
general. We note here yet another situation where, in 
the absence of the lefthandedness property, a high po- 
tential barrier might appear. For this it is convenient to 
consider a more general setup where the spatial area 2 
in Fig. [T] is filled with a medium of permittivity £2(1-1-') 
and permeability //2(w). Then, apart from a scaling fac- 
tor due to local-field correction, 9£2('^io)/[2e2(^io) + 1]^ 
|26| , the vdW potential retains the form of Eq. ([31]) , with 
the reflection coefHcicnts (15^ and p3p being modified 



as follows: \e\ — 1 (|^| — 1) is replaced by |e| — |£2| 
(I^P — lA*2p) and £ ± 1 (^ ± 1) is replaced by £ ± £2 
(^±^2). If we choose |£| > |£2| and [/i| < |^2|, the 
first terms, rather than the second terms, in the modified 
equations (|32p and (|33p dominate, but it remains that p- 
(s-)polarized waves give rise to an attractive (repulsive) 
potential. In this scenario, the appearance and the height 
of the barrier are determined by the difference between 
the absolute values of the permittivity on the two sides 
of the interface and the accordant difference between the 
absolute values of the permeability. Negative (real parts 
of the) permittivities and/or permeabilities are allowed 
but not a prerequisite. The presence of a potential bar- 
rier in this case is essentially a surface effect, the finite 
thickness of the slab and the perfect mirror playing only 
secondary roles. 

Let us return to the influence of the thickness of the 
LHM slab. As can be inferred from Eq. ([33|) . the mag- 
nitude of |Rer2_| is about 1 for d ^ (slab absent), and 
is typically determined by a e^"^ term otherwise. There- 
fore the presence of the slab is crucial for the appearance 
of a potential barrier. When the slab is very thick, the 
influence of the mirror vanishes e"^"^"^ — > [cf. Eqs. ((32|) 
and ([55)1 ]. and it is not difficult to verify that Eq. (IFll) 
reproduces the result for the resonant part of the poten- 
tial of an excited atom in front of an interface (see, e.g., 
Ref. [iy|), which, in the nonretarded limit reads as 

U{ZA) = 3-7— TTT? (37) 
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[C = [Id^op + 2|d];on/(327r£o)]. Clearly, when the slab 
is too thick, then the atom only feels the vacuum-slab 
interface. Note that a repulsive behavior of the poten- 
tial is also possible in the case of non-magnetic mate- 
rials provided that |£| < 1. It is particularly apparent 
that |J7(z^)| cannot become infinity as za —>■ 0, because 
the underlying macroscopic quantum electrodynamics is 
valid only for atom-surface distances much larger than 
the average distances between the constituents of the me- 
dia involved. 

Although potential barriers in planar magnetodielec- 
tric multilayer structures have also been found in the 
case of ground-state atoms [l^, it is worth noting that 
the barriers are generally much more significant in the 
case of excited atoms. For instance, the peaks of the po- 
tentials in the numerical examples given here are about 6 



orders of magnitude larger than those given in Ref. [l5| . 
This is because the potential in the case of ground-state 
atoms is an integral over the whole frequency range, for 
most of which the material does not exhibit strong elec- 
tric or magnetic responses. The excited atoms essentially 
allows one to pick out a frequency of choice. 



IV. SPONTANEOUS DECAY RATE REVISITED 

Whereas the (resonant) vdW potential depends on the 
real part of the Green tensor, the rate of spontaneous de- 
cay is determined by its imaginary part. From Eq. ([221) 
for the Green tensor of the absolutely nonabsorbing setup 
considered in Sec. IIII A[ it can be seen that the contri- 
butions from evanescent waves, which give rise to diver- 
gences in the region za < d, are purely real and thus do 
not contribute to the decay rate. The decay rate is thus 
expressed in terms of traveling- wave contributions which 
are finite for arbitrary za- However, baring in mind the 
findings for the vdW potential, one should carefully ex- 
amine whether these idealized results valid for nonab- 
sorbing material are an appropriate approximation to the 
reality. 

To do so, we note that from Eq. ([6]), the decay rate 
can be rewritten as 



= 1- 
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(39) 



is the decay rate of the atom in free space. Let us 
first assume an absolutely nonabsorbing LHM having 
£ — ^ — —\ again. From Eqs. (|23p and ([M)l we obtain 

Im GW = £i|j [sin(z) - ~z cos(z) - ~z^ sin(f )] , (40) 



ImG^V = 1^^^^ [sin(5) - 5cos(z)] 



(41) 



[z = 2ijJin{zA — d)/c\. Equation p8)) together with 
Eqs. ((40|) and (|4T|) . which mathematically holds for any 
atom-surface distance, is studied in Ref. [20|. It is not 
difficult to see that ImGi:^ and Im Gzz are even func- 
tions of z, and are finite at the surface. The position 
ZA = d, termed focus point in Ref. [20] , stands out in that 
the optical path for a ray of light for a round trip to the 
mirror and back is zero. For za = d, Eq. ([55]) together 
with Eqs. ([40)) and (|4ip implies complete inhibition of 
spontaneous decay (F — 0) for a dipolc moment oriented 
parallel to the surface and enhancement of spontaneous 
decay (r = 2ro) for dipole moment oriented perpendicu- 
lar to the surface [20| . 

It is worth noting that absorption can drastically 
change this behavior. Inclusion in the calculations of 



the effect of material absorption requires the calcula- 
tions to be performed on the basis of the exact (scat- 
tering part of the) Green tensor as given in Eq. (fT8|l to- 
gether with Eqs. ([Sjl- fTU]) . Numerical examples are given 
in Fig. [5l where the case of zero absorption in accordance 
with Eqs. ((38|) -(|4T]) is also shown to facilitate compari- 
son. We see that, whereas in the case of strictly zero 
absorption the decay rate as a function of the atomic po- 
sition ZA > is symmetric with respect to the position 
Za = d, any absorption destroys this symmetry. As a 
result, large enhancement of the spontaneous decay can 
be observed when the atom is near the surface, which is 
obviously due to the absorption-assisted atomic coupling 
with evanescent waves — an effect already known for or- 
dinary materials (see, e.g., Ref. [2^) which features qual- 
itatively new distance terms [cf. Eq. (1451)1 . This result is 
consistent with those reported in Ref. [22 • In particular 
for distances za < d, Eq. ([55]) together with Eqs. ([301) 
and ([1T|) can not be regarded as an acceptable approxi- 
mation to the spontaneous decay rate in the case of small 
absorption — a result which corresponds, for za < d, to 
the result found in the case of the vdW potential. 

To further elucidate the influence of the evanescent 
waves, let us examine the near-surface limit of the rate of 
spontaneous decay. By using approximations similar to 
those in Sec. lIIICl it can be shown that for za(^io/c ^ 1 
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Note that in this case only evanescent waves contribute 
to the rate. Unlike the real parts [cf. Eqs. ([5^ and ([55)) ]. 
the imaginary parts of r2„ and r2_ have the same (pos- 
itive) sign. The two polarizations therefore contribute 
constructively to the spontaneous decay rate. If the slab 
becomes sufficiently thick, Eq. ([42)) reduces, in leading 
order, to 



= 1 



D Ime 

4kTTr 



(45) 



FIG. 5: Atom-surface distance dependence of the decay rate 
of an excited atom in the setup in Fig. [T] for e = /i = —l+irj, 
d — 5c/iUio, and (a) dipole moments parallel to the surface 
and (b) dipole moments perpendicular to the surface. Curves 
1-5 refer to various degrees of absorption rj — 10~^, 10~^, 
10~*, 10~^ and 0, respectively. 



spontaneous decay, as illustrated in Fig. [S| The figure 
shows that almost complete inhibition of spontaneous de- 
cay is possible for both parallel and perpendicular orien- 
tations of the dipole moment with respect to the surface 
of a slab that is not lefthanded. 



]d!oP 



2|dJ-o|2]/(8a;3o|dio|2) ], i.e., the spon- 



[D = 3c3 

taneous rate for an atom in front of an interface. Equa- 
tion ([l5)) shows that the decay rate takes on large values 
as ZA ^0, which is a consequence of direct energy trans- 
fer from the atom to the constituents of the medium (see, 
e.g., Ref. HI). 

It is well known that lefthandedness of the slab is of 
course not a necessary condition for, say, inhibition of 



V. CONCLUDING REMARKS AND SUMMARY 

The macroscopic theory employed in this paper breaks 
down on length scales which are comparable to the size of 
the elementary building blocks of the metamaterial em- 
ployed. In order to observe the predicted effects, one thus 
has to ensure that both the atomic transition wavelength 
and the atom-surface separation are larger than this 







ZAO^lo/c 



FIG. 6: Atom-surface distance dependence of the decay rate 
of an excited atom in the setup in Fig. [T] for d = 5c/a;io, 
and (1) s — 1, fi = —1 + ilO~^, and atomic dipole moment 
perpendicular to the surface; (2) £ = —1 + ilO^'', /i — 1, and 
atomic dipole moment perpendicular to the surface; (3) e = l, 
f^ = —1 + ilO~^, and atomic dipole moment parallel to the 
surface. 



length scale. With currently available metamaterials, 
this can most readily be achieved with polar molecules 
whose rotational and vibrational transition wavelengths 
can be very large. For example, the bulk meta-material 
reported in Ref. Q exhibits negative refraction with 
Ren w —1.6 at a frequency of 1.03 THz; it consists of 
unit cells with lattice constant a « 6 x 10~^m. The fre- 
quency of the first rotational transition of a NH molecule 
(about 1 THz 28]) falls within this frequency window, 
with the associated wavelength being sufficiently large in 
comparison with the lattice constant. One could expect 
the macroscopic results to be valid down to molecule- 
surface separations of za = 10^'^m, which corresponds to 
ZAijJnk/c « 2. Smaller values of ZAOJnk/c can be reached 
by using molecules with larger transition wavelengths. 

Let us briefly comment on the strength of the force 
F{za) = —dU{zA)/dzA associated with a potential of 
the type shown in Figs. [3] and [H In this context, 
it may be convenient to express the potential and 
the force in terms of the spontaneous decay rate F. 
From Eqs. ®, (gD, and IMD it follows that U{za) = 

— 3/2fi,FoAioRe Ga;a; (^A, ^A, ti^io); where a dipole moment 
parallel to the surface has been assumed. Consider an ex- 
cited hydrogen-like atom with a transition having spon- 
taneous decay rate Fq '^ 6 x 10^ s^^ and wavelength Aio ^ 
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If the atom is at rest relative to the slab, 



the force that can counter the gravitation force and so 
levitate the atom has to be larger than mg — i x 10~^^ N. 
This can be achieved easily. For example, the peak value 
of the force corresponding to the potential barrier repre- 
sented by the dotted line in Fig.El^a) is lO'^^N. 

Consider now a situation where one has a sample of a 



dilute gas of excited atoms at some temperature. Take 
again hydrogen-like atoms with the same excited level as 
above. At a temperature of T= 10 K, the kinetic energy 
of each particle is E = 5kT/2'^2 x lO'^^ J. A LHM slab 
of thickness dujio/c = 9, permittivity and permeability 
e = /i=— l + i2x 10~^ then provides a potential barrier of 
a height of [/ = 3 x 10~^^ J. This potential barrier is high 
enough to levitate the sample of the dilute atomic gas. 
With an appropriate arrangement, it clearly can be em- 
ployed as a trapping mechanism. One may also think of 
an (evaporative) cooling apparatus. The kinetic energy is 
an averaged quantity. Some atoms among the gas sample 
can have kinetic energy exceeding the average and higher 
than the potential barrier. These can overcome the bar- 
rier and get adsorbed by the walls, in so doing lowering 
the average kinetic energy of the remaining (a little more 
diluted) gas sample. It should be pointed out that in the 
examples meta-materials with very small absorption are 
assumed, which are still beyond the reach of today's ex- 
perimental techniques. Furthermore, since excited atoms 
are considered, the effects predicted hold only for time 
periods short compared to the atomic lifetime. A more 
refined theory should include the broadening and shifting 
of the atomic transition frequency. 

In summary, we have studied the resonant vdW po- 
tential experienced by an excited two-level atom in front 
of a meta-material slab backed by a perfect mirror, with 
emphasis on a LHM slab with unity negative real part 
of both the permittivity and the permeability, and an 
atom located in vacuum. We have shown that material 
absorption has to be taken into account in order to avoid 
divergent values of the potential when the atom-surface 
separation becomes smaller than the slab thickness. Only 
for values of the atom-surface separation which are larger 
than the slab thickness, disregarding of absorption leads 
to reasonable results provided that absorption is suffi- 
ciently small. 

Further, we have shown that the setup can be used 
to form potential barriers for excited atoms - barriers 
that are typically of several orders of magnitude higher 
than those observed for ground-state atoms. In particu- 
lar, they are significant enough to suggest potential use in 
quantum levitation, trapping, or cooling of atoms. The 
barriers are results of competing contributions from s- 
and p-polarized waves, thus are peculiar to atoms with 
dipole moments parallel to the surface. For random 
atomic dipole moment orientations, a corresponding av- 
eraging has to be applied which may lead to a reduction 
of the barrier height. 

Finally, we have revisited the effect of the LHM slab 
on the spontaneous decay of the excited atom. We have 
shown that already an arbitrarily small but finite amount 
of material absorption can drastically change the decay 
compared to that predicted when absorption is ignored. 
In particular near the surface, absorption gives rise to 
a noticeable enhancement of the decay. As a prospec- 
tive task, the theory could be expanded to other systems 
such as spherically or cylindrically symmetric ones, am- 
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plifying instead of absorbing media could be assumed, 
and the vdW interaction between excited atoms near a 
macroscopic body of meta- material could be investigated. 
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